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Metrological investigations of the quantum Hall effect (QHE) completed by transport measure- 
ments at low magnetic field are carried out in a- few- /xm- wide Hall bars made of monolayer (ML) or 
bilayer (BL) exfoliated graphene transferred on Si/SiOj substrate. From the charge carrier density 
dependence of the conductivity and from the measurement of the weak localization quantum correc- 
tion at low magnetic field, we deduce that transport properties in the devices are mainly governed 
by the Coulomb interaction of carriers with a large concentration of charged impurities. In the QHE 
regime, at high magnetic field and low temperature (T < 1.3 K), the Hall resistance is measured by 
comparison with a GaAs based quantum resistance standard using a cryogenic current comparator. 
In the low dissipation limit, it is found quantized within 5 parts in 10^ (one standard deviation, 
la) at the expected rational fractions of the von Klitzing constant, respectively Rk/2 and Rk/'I 
in the ML and BL devices. These results constitute the most accurate QHE quantization tests in 
exfoliated graphene to date and the first in BL graphene. It turns out that a main limitation to the 
quantization accuracy, which is found well above the 10~^ accuracy usually achieved in GaAs, is 
the low value of the QHE breakdown current being no more than 1 /lA. The current dependence of 
the longitudinal conductivity investigated in the BL Hall bar shows that dissipation occurs through 
quasi-elastic inter-Landau level scattering, assisted by large local electric fields. We propose that 
charged impurities are responsible for an enhancement of such inter-Landau level transition rate 
and cause small breakdown currents. 

PACS numbers: 73.43.-f, 72.80.Vp, 06.20.-f 
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I. INTRODUCTION 

The discovery of the quantum Hall effect (QHE) in 
1980i has revolutionized resistance metrology by estab- 
lishing a universal quantum resistance standard at ra- 
tional fractions of the von Klitzing constant Rk = h/e'^ 
where e is the electron charge and h Planck's constant. 
Although the QHE was first observed in Si-MOSFETs, 
the cleaner two-dimensional electron gas (2DEG) made 
by epitaxial growth of GaAs/AlGaAs heterostructure 
provided more practical quantum resistance standards. 
They give accurate and reproducible representations of 
i?K within an uncertainty below one part in 10^ when 
operated at low temperature (T = 1.5 K) and high mag- 
netic induction {B = 10 T)^i^. Following observation of 
the QHE in graphene^i^ with a sequence of Hall resis- 
tance plateaus at i?H = ±-RK/(4(n + 1/2)) (with n an 
integer ^ 0) that survive even at room temperature^, an 
application to resistance metrology was considered^. The 
peculiar QHE originates from the honeycomb lattice of 
carbon atoms in which charge carriers at low energy be- 
have like chiral massless relativistic fermions with Berry's 
phase 7r~. Under magnetic field, the density of states be- 
comes quantized in Landau levels (LLs) with a AeB/h 
degeneracy (valley and spin) that occurs at energiesS, 
zbwF yJ^hneB / c. The robustness of the QHE on the first 
plateau comes from the energy spacing 36-\/B[T] meV 
between the first two LLs being larger than in GaAs 
(1.7i?[T] meV). In bilayer graphene, which consists of 



two graphitic monolayers with Bernal stacking, the dis- 
persion relation becomes parabolic and carriers behave 
like chiral massive (m = 0.033 x mg with mg the elec- 
tron mass)^ Dirac fermions with Berry's phase 27r — . 
This leads to an integer QHE"' with resistance plateaus 
at i?H = ±^^K/(4(n + 1)), with n an integer ^ 0. The 
energy gap between LLs occuring at^^ ±hjJc^J n{n — 1) 
[lJc = eBjm is the cyclotron pulsation) is smaller than 
in single graphene layer, especially at low magnetic field, 
but is larger than in GaAs systems. Larger energy gaps 
give much hope that a more practical resistance standard 
operating at a lower magnetic field or a higher tempera- 
ture could be developed in both graphene systems. In the 
short term, comparison of the Hall resistance in graphene 
systems and in GaAs would constitute a stringent test of 
the QHE universality. This would support ongoing efforts 
to make an historic evolution towards a Systeme Interna- 
tional of units directly linked to fundamental constants 
of physicsiS,. More generally, the metrological approach 
can supplement the understanding of physics to the lim- 
its of instrumentation. Lastly, meeting the very demand- 
ing metrological requirements for the QHE application in 
graphene (quality of electrical contacts, control of elec- 
tronic properties such as mobility and density over large 
mm-size scale) further enhances the severeness of the 
benchmark test offered by the QHE for the quality of any 
two-dimensional material, and makes it very significant 
and useful for the development of industrial applications 
such as microelectronics. The metrological investigation 
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has started shortly after the discovery of the QHE in 
graphene. Previously, the Hall resistance i?H was demon- 
strated to agree with i?K/2 on the plateau corresponding 
to Landau level filling factor v — nsh/{eB) = 2 in exfo- 
liated monolayer graphene within a relative uncertainty 
of 15 parts in 10^ (one standard deviation, la), probably 
limited by the high resistance of contacts (> 1 kO)i^. 
More recently, Tzalenchuk and co-workers have reported 
an agreement within an uncertainty as low as 9 parts in 
10^^ (Icr) in a large sample (160 x 35 //m^) made of epi- 
taxial monolayer graphene grown on the Si-terminated 
face of silicon carbide (SiC), with a mobility of about 
7500 cmVV/s when placed at B = 14 T and T = 0.3 K 
— ii^. Achieving the QHE quantization in graphene with 
similar uncertainty at a few teslas magnetic induction 
and higher temperature, which is required to develop a 
quantum resistance standard challenging the GaAs ones, 
is still a critical issue. 

In this paper, we report on the accurate investiga- 
tion of the QHE quantization in exfoliated graphene on 
Si/Si02 substrate, for the first time in a bilayer (BL) and 
also in a monolayer (ML). Measurements were performed 
with a Cryogenic Current Comparator (CCC)-based re- 
sistance bridge. The objective was to determine the lim- 
itations to the quantized Hall resistance accuracy previ- 
ously experienced in exfoliated graphene, which turned 
out to be the reference technique enabling to unveil most 
of chiral Dirac fermions electronic transport properties. 
The understanding of these limitations could even be use- 
ful to overcome likely obstacles in the development of 
quantum resistance standards with higher performances 
in graphene grown either on SiC or by Chemical Vapor 
Deposition (CVD). 

The paper is organized as follows. In Section II, we 
report on electronic transport properties of graphene in- 
vestigated by means of conductivity measurements at low 
magnetic field. In both the ML and BL samples, the 
analysis of the conductivity dependence on charge car- 
rier density shows that carriers are mainly scattered by a 
large concentration of charged impurities located about 
1 nm close to the graphene flakes. The major impact of 
charged impurities responsible for strong spatial fluctu- 
ations of the carrier density which survive at finite den- 
sity is also confirmed by measurements of quantum cor- 
rections to conductance (weak localization and universal 
conductance fiuctuations) in the BL sample. Section HI 
reports on quantization tests performed by means of com- 
paring the QHE in GaAs and in graphene systems. For 
monolayer and bilayer graphene, the Hall resistance of 
the first plateau (Landau levels are spin and valley de- 
generated) in the zero dissipation limit is found quantized 
within 5 parts in lO"^ (Itr) to i?K/2 and -Rk/4 respec- 
tively. This is the first QHE quantization test in bilayer 
graphene. One main limitation to accuracy is the low 
value of the QHE breakdown current limited to about 
~ 1 fiA. In section IV, we show that the mechanism of 
dissipation (or backscattering) in the BL sample, which 
ends up in the QHE breakdown is based on quasi-elastic 



inter-Landau level scattering (QUILLS) assisted by large 
local electric fields. This leads to discussing the role of 
charged impurities in enhancing inter-Landau level tran- 
sitions. 
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FIG. 1. (Color online) Optical images of the BL sample a) and 
of the ML sample b) with contacts resistance values indicated 
below. 



II. ELECTRONIC TRANSPORT PROPERTIES 
AT LOW MAGNETIC FIELD 

Measurements were carried out on 15 x 2 /im^ and 
26 X 4.6 /im^ Hall bars based on monolayer graphene 
(ML) and bilayer graphene (BL) respectively, which have 
been mechanically exfoliated from natural graphite (see 
Fig. [T]). Flakes were transferred on top of highly doped 
silicon substrates covered by 90 nm (resp. 500 nm in 
BL) of thermally grown Si02 used for backgating. Sam- 
ples were finally covered with a 300 nm-thick polymethyl 
methacrylate (PMMA) resist layer. Transport properties 
were explored by four-terminal resistance measurements 
defined by i?ij,ki = {Vk — Vi)/Ii^i, where Vi is the volt- 
age potential at terminal i and is the current flowing 
between terminals i and j . 

A. Influence of charged impurity scattering on 
conductivity 

In both samples, the four-terminal conductivity a — 
l/p ~ 1/-Rij,ki X {W is the sample channel width, dki 
the distance between terminals k and 1), deduced from 
^^06. 23 a-iid -Ri8,24 measurements in the ML and BL sam- 
ples respectively, was analyzed at zero magnetic field as 
a function of the gate voltage Vg ■ It shows a typical min- 
imum that occurs at Vcinin (see Fig. [5^). At this value 
the carrier density defined as Us — Cg(Vg — VGmin)/e 
(with Cc/e = 2.40 X 10" cm'VV for ML and Cc/e = 
4.31 X lO^'^ cm^^/V for BL) is zero on spatial aver- 
age, n — — CcVGniin/e is the carrier density induced in 
the graphene by surrounding charged impurities. While 
annealing the samples under vacuum at a temperature 
of about 400 K, the conductivity dip becomes sharper 
and its position Vcmin shifts near zero indicating an in- 
crease in the carrier mobility fj, and a decrease of At 
T = 1.3 K for the ML sample (resp. T = 0.35 K for 
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the BL sample) and at carrier density away from the re- 
gion of the minimum conductivity, (t{Vq) is quite hnear 
for ML with no proof of subhnearity in the considered 
density range (< 2.5 x 10"'^^ cm~^), and is shghtly super- 
hnear for BL. These features indicate that the long-range 
Coulomb potential induced by charged impurities consti- 
tute the dominant source of scattering in the considered 
samples^^'^^. 

Conductivity for the ML sample, except near the min- 
imum {—n ± 5 X 10^^ cm~^), is well fitted by the theo- 
retical model based on Boltzmann transport theoryi^iiS 
a^CQVo/e) = a{ns — n) — G(rs,d)^^ valid for elec- 
trons (ris > n*) and for holes (ns < —n*) where rii 
is the density of charged impurities at an average dis- 
tance d from the conductor (in the silicon substrate or 
in the PMMA). n* is a residual density corresponding 
to the density of electron and hole puddles into which 
the system breaks at low density because of the inho- 
mogeneous density profile created by Coulomb impuri- 
ties, describes the full dielectric environment of the 
sample that screens Coulomb interactions. Considering 
two semi-infinite media made of Si02 and PMMA on 
top of the device with dielectric constants esi02 = 3.9 
and epMMA = 4.5 respectively, = 2e^/(47reo(epMMA + 
esi02)?iWF) = 0.47 (with ?;f = 1.1 x 10^ ms'^ ^) and 
G{rs — 0.47, d = 0) = 28.2 in the random phase ap- 
proximation. It appears that the dielectric constant of 
PMMA higher than air or vacuum screens more efficiently 
the coulombic potential of charged impurities. Note that 
since G{rs,d) is only weakly dependenti^ on d, the ap- 
proximated value G{rs,d = 0) is valid while the elec- 
tron/hole asymmetry in the conductivity curve remains 
weak, as observed, and thus is not considered. The mean 
impurity density (electron/hole average) deduced from 
the adjustment is rii w 1.9 x 10^^ cm~^. At low density, 
assuming this value of n; and a finite value of d in the 
range lower than 2 nm, the Boltzmann transport theoryi^ 
correctly explains (within a factor of 2, see ref<^) the ex- 
perimental values of the conductivity minimum ao, of 
the plateau width minimum conductivity n* and of the 
minimum position —n. The size ^ and the density n* 
of electron/hole puddles near the charge neutrality point 
(CNP) can be calculated from^^ ^ ^ l/{r^^^) ^ 32 nm 
and n* = cr^n-JGir^ = 0.47, d = 0)^ = 2.7 x 10" cm^^ 

(with CTo = 4^) respectively. One deduces that each 
puddle contains about 9 elementary charges in average. 
The theoretical model can also explain the conductivity 
curve asymmetry which corresponds to a constant mobil- 
ity {fj, = a/{nse)) higher for holes (4050 cm^/V/s) than 
for electrons (3400 cm^/V/s) by a typical < -t-5 A-size 
shift of the distance d of charged impurities from the 
graphene layer under the electric field effect produced 
by the back gate voltage, assuming unequal numbers of 
random positively and negatively charged impurities!^. 
Such a similar electron/hole asymmetry has already been 
observed in dirty samples^i^. On the other hand this 
asymmetry cannot be explained by the theory for attrac- 



tive vs. repulsive scattering of massless Dirac fermions 
by Coulomb impurities^.- predicting a higher mobility for 
electrons for a negative value of n. Neither can a lo- 
cal doping due to the presence of metallic contacts on 
graphene account for it since they are non-invasive in the 
studied samples^^ and would have induced sublinearity of 
the a[CoVole) curve. 

In the BL sample, the conductivity can also be well fit- 
ted by a similar transport theory^! based on Coulomb in- 
teractions with charged impurities a{CQVQl e) = a{ns — 

") « ITX^fl + iBVWKd + Qtf)] valid for clec- 
trons {us > n*) and holes (ng < —n*) with qtf^^ = 
47reo(epMMA + esi02)fi^/8TOe^ = 0.6 nm the Thomas- 
Fermi screening length. This theory therefore includes 
the superlinearity which results in a mobility depending 
on carrier density. Not only does the conductivity adjust- 
ment give ni w 2 X 10^^ cm~^ but also d w 1 nm. These 
values are comparable to values extracted from the ML 
conductivity curve. Taking into account that the two 
samples have been made using the same technological 
processes and the same substrates (except for the Si/Si02 
thickness) , this agreement strongly supports our descrip- 
tion of conductivity using the Boltzmann transport the- 
ory based on long-range Coulomb scatterers. It appears 
that this model correctly predicts, except for the mini- 
mum position {~n), the experimental values of CTq and of 
the plateau width minimum conductivity n* . From the 
values of = 11 nm calculated with the specific model 
developed for BL^^ and the value of n* = '^ojg^^i = 
2.2 X 10^^ cm~^ (with ctq — 5-5^), one deduces that each 
puddle contains about 8 elementary charges. An elec- 
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FIG. 2. (Color online) a) Conductivity as a function of carrier 
density controlled by the backgate voltage for the BL (green) 
and ML (blue) samples, b) Magneto-conductivity in the BL 
sample at T = 0.35 K at carrier densities in the range rig = 
-2 X 10^^ cm"^ ± Ans/2 with An^ = 3.3 x lO" cm"^. c) 
Magneto-conductivity after averaging on carrier density and 
theoretical adjustments at T = 0.35 K and T — 1.5 K. 
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tron/hole asymmetry of the conductivity is also observed. 
But contrary to the ML sample, at rtg = 2 x 10^^ cm~^ 
the electron mobility (2307 cm^/V/s) is higher than the 
hole mobility (1940 cm^/V/s) by about 16%. It can again 
be explained by the shift of the mean distance d between 
the charged impurities and the graphene layer by a few 
A under the electric field produced by the voltage on 
the backgate, but with impurities in excess with a sign 
opposite to the ML sample casein. We note that the 
same amount of charged impurities leads to a lower car- 
rier mobility in the BL sample than in the ML sample, 
confirming that long-range Coulomb scattering is a very 
efficient mechanism to spoil mobility in bilayer graphene. 
Beyond providing a very efficient source of scattering, 
charged impurities give rise to strong spatial fluctuations 
of the charge carrier density with a correlation length 
^ of 32 nm and 11 nm for ML and BL samples respec- 
tively. These fluctuations leading to electron/hole pud- 
dles landscape near the CNP are also known to persist 
at the higher carrier densities (rtg ~ 1 — 2 x 10^^ cm^^) 
where the QHE has been investigated32i^. Also, a more 
macroscopic inhomogeneity of the carrier density at a 
/im-size scale with a typical amplitude of a few lO^^cm"^ 
has been observed. In particular, it manifests itself in 
the BL sample through spatial variation of —n, ao and 
of da/dVQ slopes. These quantities depend for instance 
on the conductor area probed in different configurations 
(e.g. i?i8,24 and i?i7,34). These carrier density fluctua- 
tions can be explained by spatial variations of n; by a 
few 10^^ cm~^ and of d by a few A, which is also the 
typical height of graphene flake ripples. Therefore the 
samples are far from being homogeneous compared to 
GaAs based 2DEG commonly used for quantum resis- 
tance standards, where less than 10^" cm^^ variation of 
Us can be achieved. 



In both samples, the diffusion coefficient D and trans- 
port mean free path ^tr can be determined from conduc- 
tivity measurements at low temperature and carrier den- 
sities where the QHE was investigated. In the ML sam- 
ple, at r = 1.3 K, and rig = 6.4 x 10^^ cm~^ (electrons), 
corresponding to a Fermi energy of Ep = hv-p^/n^ = 
102 meV, D is calculated using the Einstein relation 
D = a{ns)y/TThvY / {2e^ = 2.0x10"^ m^s"^ and then 
Ztr = 2D/vf = 36 nm. In the BL sample, at T = 0.35 K, 
and rig = —2 x 10^^ cm~^ (holes), i.e. Ep = hy^n^/m — 
72 meV, D = a{ns)nh'^ /{2e^m,) = 1.5 x 10"^ ^2^-1^ 
Itr = 2D/vp{ns) = 34 nm, with up (rig) — hy/n^/m — 
8.8 X 10^ ms~^. These values confirm that electronic 
transport is diffusive with similar amount of disorder 
in both samples: kplu- — 5.2 for ML and 8.5 for BL. 
Comparatively, in cleaner GaAs/AlGaAs Hall bars used 
as quantum resistance standards, with typical mobili- 
ties 280000 cm^/V/s and density 5.2 x 10" cm'^, /j, = 
3.4 fj,m is 100 times higher and kpltr ~ 600. 



B. Quantum corrections to conductivity 

In the BL sample, where the QHE has been 
more extensively studied, quantum interference correc- 
tions to conductivity, both weak-localization correction 
(WL) and reproducible mesoscopic conductance fluctu- 
ations (CF), were investigated by performing magneto- 
conductivity measurements at low temperature. Actu- 
ally, both in monolayer and bilayer graphene the am- 
plitude of these corrections is not only ruled not by 
inelastic scattering like in any other diffusive metal, 
but also by elastic scattering mechanisms affecting the 
valley symmetry (intravalley scattering and/or trigo- 
nal warping of the conical band structure, interval- 
ley scattering). This is a consequence of the direct 
manifestation of chirality property in quantum interfer- 
ence effects. For instance, interferences between time- 
reversal symmetric diffusive electron trajectories lead to 
weak-localizatioiii^orrections to conductivity in bilayer 
graphene because of the charge carrier wave function 2tt 
Berry's phase while weak-antilocalization^is expected in 
monolayer due to tt Berry's phase. 

Applying a magnetic field breaks the system time re- 
versal symmetry and suppresses the weak-localization 
corrections. This gives rise to a well-known magnetoresis- 
tance. Four-terminal magnetoresistance measurements 
R{B) were carried out at temperatures T = 0.35 K and 
T = 1.5 K, at low magnetic field using a standard AC 
low-frequency (13 Hz) lock-in technique, and a low-noise 
preamplifier. The measurement current is / = 30 nA, 
thus the effective temperature of carriers assessed by 
Toff = eRI/k-Q = 0.52 K, where R is the resistance 
per square, is slightly higher than the base temperature 
0.35 K. Fig. [2}d reports a set of magnetoconductivity 
curves recorded at densities around ns = — 2x10^^ cm^^ 
over a total range Arts = 3.3 x 10^^ cm~^ . They all 
display a characteristic dip at zero field, signature of 
the expected weak-localization, the amplitude of which 
barely exceeds reproducible fluctuations (CF) which are 
analyzed below. To make WL conductivity deep stand 
out from fluctuations, magnetoconductivity curves were 
averaged over the full density range where the diffusion 
coefficient D does not vary by more than 10%. The aver- 
aged curve (see Fig.[2t) is then adjusted by the appropri- 
ate weak-localization theory^^^, A(j{B) — a{B) ~ (t(0) — 

|i[f-(_:^)„f( Z^-;^ ) + 2f( ^_, )]. Here 

F{z) — ln{z) + -0(1/2 + z^^), ■ip{x) is the digamma func- 
tion, Tg^ = 4eDB/h. T^\JyL ~ ^/^IwL is the phase 
breaking rate. t.~^ = D/Lf is the intervalley scattering 
rate lifting the valley degeneracy of electronic states and 
which is caused by short-range defects with maximum 
size of the order of the lattice spacing, t^^ — D/L^ — 
2t^^ + t^^ is an intravalley scattering rate, t^^ is the 
intravalley chirality breaking rate caused by surface rip- 
ples, dislocations and atomically sharp defects, i.e. short- 
range defects. is the intravalley p — —p symmetry 
breaking rate (where p = hkp, kp is the carrier momen- 
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turn at the Fermi level) caused by the anisotropy of the 
Fermi surface in k space, i.e. the trigonal warping. In 
bilayer graphene, assuming a quadratic Hamiltonian, it 
is expected that = t^.^ where Ttr = hr/vp is the 
transport timc"^^. 

The adjustments of data at T = 0.35 K and T = 1.5 K 
give the phase coherence length L^y^i^iT = 0.35 K) = 
0.47 ^m and L$wl(T = 1.5 K) = 0.42 /xm, the in- 
tervalley scattering length Li(r — 0.35 K) = 0.51 /im 
and Li(T — 1.5 K) = 0.47 /im, the intravalley scat- 
tering length L^{T = 0.35 K) and L^{T = 1.5 K) 
< 0.03 /im. The extracted values are very similar to 
those measured in bilayer graphene and reported in the 
literature^. -L$wl below the sample size indicates that 
electronic transport is not fully quantum coherent. It 
appears that L$wl ~ L\ and Li 3> i*. The WL is 
made observable due to significant intervalley scattering, 
though much less than intravalley processes. It also ap- 
pears that ~ Ztr- The fact that = 2t~^ -|- ~ 
T^^^ means that t^^ is small, since it is expected that 
= Tt~^. Finally, the fact that Li ^ L^. ~ Itr, to- 
gether with <C demonstrate that short-range 
scattering is not dominant. Moreover, L$wl appears 
quasi constant between T = 0.35 K and T = 1.5 K 
around a value that is far below the typical size of the 
sample. Such saturation of L$wl at low temperature, 
well below the particle-particle interaction length (Lhh — 

V'^[ 2.e^MW(2e^) T5^ = 1^^ T = 0.35 K)M, has 
already been observed in graphene samples^ near the 
CNP. It could be a feature of transport by percolation 
through electron/hole puddles"^^ persisting at finite den- 
sity (typically rig = —2 x 10^^ cm^^ < n*) in the very 
inhomogeneous BL sample. These results confirm the 
conclusion drawn from the analysis of the conductivity 
curves a{CQ,VQ,/e) that long-range Coulomb scattering 
by charged impurities trapped in the silicon substrate or 
in the PMMA top-layer of the graphene-based sample is 
dominant. 

Conductance fluctuations were measured by varying 
the magnetic induction over a ±1 T magnetic field range 
with a measurement current of 50 nA at T = 0.35 K. 
The standard deviation is found to be 5Gb = 0.021e^//i. 
In graphene, CF resulting from interference of phase- 
coherent chiral carrier diffusive paths are also expected 
to depend on elastic scattering. In the BL sample, 
since Z/4.WL ~ L^, one expects the amplitude of CF 
to be properly described by the theory of well-known 
universal conductance fluctuations (UCF) for diffusive 
metals^i^. Precisely, in the case of a two-dimensional 
conductor, at a magnetic field larger than the typical 
magnetic field of WL magnetoresistance, it is given'^^ by 
5G = 0.862-^5^^ ^jjg^g ^ ^ 

the leng th and th e width of the conductor measured, 
Lt — yJhD /k-QT is the thermal length. Assuming the 
values W = 4.6 ^im, L ^ 8 ^m, T^s = eRI/kB = 0.87 K, 
one finds L^^^^ = 0.36 /im < = 0-47 /im which 

results in SG = O.Q21e^/h. The good agreement of 



the experimental magnitude of CF with the theoreti- 
cal value of the UCF in diffusive metals confirms that 
Li ^ L$ ~ 0.5 /im, and since Ztr = 34 nm <C Li, that 
long-range scattering is dominant. On the other hand, it 
shows that conductance fluctuations as a function of the 
magnetic field are not sensitive to the observed carrier 
density inhomogeneity or presence of electron and hole 
puddles. 

The analysis of transport at low-magnetic field shows 
that the dominant mechanism of scattering in our sam- 
ples is Coulomb interaction with a large concentration of 
charged impurities closely surrounding graphene fiakes 
(in the silicon substrate and in the PMMA top layer cov- 
ering the devices). Beyond to drastically reducing the 
carrier mobility, they are responsible for strong spatial 
fluctuations of the carrier density that might stay bipo- 
lar even at finite density (a few 10^^ cm~^). 

III. HALL RESISTANCE QUANTIZATION 
TESTS OF THE QUANTUM HALL EFFECT 
REGIME 
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FIG. 3. a) Hall resistances (at magnetic inductions B = 
2, 5, 7, 9, 13, 18.5 T) and longitudinal resistance (at 
B = 18.5 T) in the BL sample at T = 0.35 K. b) i?xx = 
fii8,24 X ^ (blue), _Rxx = -R28,34 X ^ (magenta) for cur- 
rents I = 0.5, 1, 3 fiA. c) i?xx ~ i?i8,24 X ^ as a function 
of Us around u = —4 at T = 0.35 K and I = 0.5 fiA (blue); 
T = 0.35 K and / = 1 /iA (deep blue); T = 1.5 K and 
I = 0.5 /iA (orange); T = 1.5 K and / = 1 /lA (red). Ver- 
tical dashed lines underline reproducible fluctuations. Inset: 
/n(i?xx) as a function of the carrier density, d) Three-terminal 
resistance of contacts as a function of Hs. 



In the QHE regime, all measurements were performed 
using direct current (DC) measurements techniques. Re- 
sistances values reported in the following are the average 
of values measured for both current directions. About 
notations, i?xx is a longitudinal resistance value normal- 
ized to a square, for example -Rxx = ^ij,ki ^ ^ if the 
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longitudinal resistance is measured between terminals k 
and 1. Fig. [3^ shows Hall and longitudinal resistance as 
a function of carrier density for the BL sample. Mea- 
surements clearly reveal v = ±4 Hall plateaus becoming 
well defined at the highest magnetic inductions, v = ±8 
plateaus are barely visible. We note that the energy gap 
between the lowest LLs (n = 0, 1 and n = 2) is 92 meV 
(1068 K equivalent temperature) at B = 18.5 T, thus 
3050 times the thermal energy at T = 0.35 K. The longi- 
tudinal resistance i?xx reported at B = 18.5 T exhibits a 
central peak corresponding to the degenerate n = and 
n = 1 LLs and minima occuring simultaneously with Hall 
plateaus. We only investigated the physics of the v = —4 
plateau for holes, which is more flat and characterized 
by a drop to zero of i?xx- Dissipation level in the 2DEG 
and quality of contacts are essential quantization crite- 
ria of the QHE, as demonstrated by several experimental 
works'^^ as well as the Landauer-Biittiker theory^. The 
quantization is indeed directly related to the absence of 
dissipation {i.e. of backscattering) , the rate of which can 
be determined by the measurement of i?xx- Fig.[3]D shows 
the behavior of i?xx with hole density on the ly — —4 
plateau in the BL sample for several current values in- 
creasing from 0.5 /LtA to 5 ^A. The i?xx plateau shrinks 
and simultaneously the i?xx minimum increases. Fig. [3]d 
also shows that position and magnitude of i?xx minima 
depend on the sample region measured. Position varia- 
tion can be attributed to carrier density fluctuations with 
a magnitude of a few 10^"'^ cm~^ caused by charged im- 
purities, as already mentioned in Section H. In addition, 
the magnitude variation illustrates that the ignition of 
QHE breakdown is a very spatially inhomogeneous phe- 
nomenon. Fig. [5}; shows that the temperature effect on 
i?xx between 0.35 K and 1.5 K is smaller than the cur- 
rent effect between 0.5 fiA and 1 /nA. It also shows that 
i?xx has reproducible fluctuations as a function of with 
a similar pattern at the two different temperatures and 
currents. We will later discuss the origin of these fluc- 
tuations, particularly visible a bit away from the mini- 
mum because of a better signal to noise ratio. Averaging 
fluctuations (and noise) of i?xx around speciflc density 
values gives typical and relevant mean values of the lon- 
gitudinal resistance .Rxx- At T = 0.35 K and / = 0.5 /lA, 
-Rxx = -Ri8,24 X ^ is {2± 14) m.n and (62 ± 9) mil at 
Us = -1.88 X 10^2 and = -2.01 x lO^^ ^^-2 

respectively. These resistance values are to be compared 
with 100 ^Q, the typical value of i?xx which ensures a 
10~^ i?H accuracy in usual GaAs-based quantum stan- 
dards (LEP514ii). Contact quality was determined by 
performing three-terminal measurements of resistance. 
RsT (e- 9- ^ij,ii) g;ives the resistance value of the contact 
Rc (e. g. i) combined with a i?xx (e. g. i?kj,ii) con- 
tribution. For a good contact, the drop to a negligible 
value of i?xx (^ 1 ^) in the dissipation-less state leads 
to a flat minimum of the resistance R^t giving an up- 
per bound of Rc . As observed in Fig. |3Ji, for the good 
Ti/Au contacts of the BL sample, Rc values deduced 
from i?3T minima can be as low as 10 fl (see Fig. [TJi). 



Note that i?3T minima occur at slightly different rig val- 
ues due to the carrier density spatial inhomogeneity. i?3T 
for contact 3 does not exhibit such a flat minimum with 
a value higher than 428 il. The highest resistance value 
was found equal to 5.9 kil for contact 8. These anoma- 
lous behaviors can be explained by a large fluctuation of 
Us in the voltage arm thin channel (2 /jm) or even by a 
partial breaking of the constriction probably caused by 
the sample cooling down too fast. The complete breaking 
can account for the infinite resistance observed for some 
other contacts. Contacts 3 and 8 were used as current 
contacts, rather than voltage, for the Hall resistance pre- 
cision measurements. It was indeed demonstrated^ that 
a very resistive detecting voltage contact can lead to a de- 
viation from quantization notably because being unable 
to restore the equilibrium of the edge state population^. 
Although we used Pd instead of Ti/Au to make contact 
to graphene, similar observations are reported in the ML 
sample. The five contacts used to perform measurements 
have low resistance values ranging from 15 $7 to 260 
(see Fig. [It). 
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FIG. 4. a) Relative Hall discrepancy ARh/Rh as afunction 
of Tis (and v in upper-scale) at four currents, b) ARu/Rh as 
a function of J at ris = —2.01 x 10^^ cm~^ (blue) and at 
ris = —2.1 X 10^^ cm~^ (red), c) i?xx = -Ris 24 x ^ as a func- 

^ ' ' ' "24 

(filled blue square) and at 
(filled red circle), iixx = i?28,34 x ^ 
as a function of I at Ub — —2.01 x 10^^ cm~^ (unfilled green 
square), d) _Rxx = i?i8,24 x ^ as a function of / for filling 
factors = —4.2 (green) and v — —4.8 (blue) at tempera- 
tures T = 1.5 K and T = 0.35 K. Error bars correspond to 
uncertainties given within one standard deviation, la. 



We then performed accurate measurements of i?H in 
terms of Rk using a resistance bridge equipped with a 
SQUID based cryogenic current comparator. In practice, 
the Hall resistance is compared to a well-known 100 51 
wire resistor calibrated in terms of a GaAs based quan- 
tum resistance standard (LEP514). In the BL sample, 
Fig. 13^ reports the relative deviation of i?H = ^38,24 
from its nominal value l^Ru/Ru = i?H/(-RK/4)-l as a 
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function of Us- All uncertainties are given within one 
standard deviation (la). Let us note that the resis- 
tance measured not only includes a pure transverse re- 
sistance but also a longitudinal resistance contribution, 
because the line between voltage terminals is not per- 
pendicular to the one between current terminals. Mea- 
surements clearly show a flat resistance plateau within 3 
parts in 10^ over a 2 x 10^^ cm~^ carrier density range 
when measured with a current below 1 ^A. At the low- 
est measurement current / = 0.5 /iA, deviations from 
quantization at highest carrier density agree with the ex- 
pected shape of the Hall plateau (decrease of resistance 
on plateau edges). The shape evolution at higher cur- 
rents is attributed to a i?xx contribution which adds to 
the transverse resistance and increases with the current. 
This coupling between i?H and i?xx, which always ex- 
ists to some extend in GaAs based quantum resistance 
standard''^, will be later discussed in more details. The 
flatness appears worse at / = 2 /iA, as expected with re- 
gards to the large increase of i?xx, fluctuating with carrier 
density as previously discussed. Fig. [Jb confirms that 
deviations from quantization start to drastically increase 
from 1 — 2 fiA at rig — —2.1 x 10"'^^ cm~^ and from 
/ = 3 /iA at Hs = —2.01 X 10"^^ cm~^. As demonstrated 
in Fig. I3J:, the increase of deviation due to current is ac- 
companied by a large increase of i?xx at both densities. 
The weighted mean value of ARu/Ryl values measured 
at currents below these critical currents leads to small 
deviations of (0.57 ± 3.1) x 10"^ and (3.0 ± 3.2) x 10"^ 
at Us = -2.01 X 10" cm-2 and = -2.1 x 10" cm"^ 
respectively. Since i?xx is the relevant parameter of quan- 
tization, ARn/ Rusls a function of i?xx is then reported in 
Fig.[5]from data of|3}D and|3t for the two carrier densities. 
Although dissipation is inhomogeneous in the sample, as 
the very different values of i?xx measured using voltage 
terminal-pairs (2,4) and (3,4) at Ug — —2.01 x 10^^ cm~^ 
express again, all deviations scale quite linearly with i?xx, 
indicating a common coupling mechanism between Hall 
and longitudinal resistances. This linear relationship, 
which is usually observed in GaAs based quantum re- 
sistance standards, is generally explained in terms of an 
effective misalignment of Hall probes, either due to a 
lack of carrier density homogeneit y ^"'''^^ in the sample 
or to current flow chiral nature in finite width voltage 
terminals^^. In a good QHR sample, one usually finds 
AR}i/ R}i=aRxx/ Ru with a ~ 0.1 — 1. In our case volt- 
age terminals are really misaligned, which should lead 
to a unity coupling factor. But from slopes we deduce 
a values in the range 10~^ to 10^*, depending where 
i?xx is measured. This means that, due to inhomogene- 
ity, i?xx values are not quantitative measurements of the 
dissipation level between Hall probes 2 and 4 when the 
current flows between terminals 3 and 8. Nevertheless, 
the values as a whole give a qualitative representation 
of the dissipation current behavior in the sample. It is 
therefore justified to extrapolate Ai?H/i?H in the dissipa- 
tionless limit (i?xx= 0) at which the perfect quantization 
is expected. In this limit, at Us = —2.01 x 10^^ cm^^. 



we find Ai?H/i?H(-Rxx = 0) = (-6.62 ± 3.0) x 10"^ 
and Ai?H/i?H(i?xx = 0) = (-2.43 ± 3.7) x 10"^ us- 
ing i?xx measurement with voltage terminal pairs (2,4) 
and (3,4) respectively. Agreement of these two values 
within the measurement uncertainty corroborates our 
extrapolation protocol. At rig = —2.1 x 10^^ cm^^, 
Ai?H/i?H(i?xx = 0) = (-0.94±3.78) x lO"'^, thus the Hall 
resistance stays quantized within the measurement un- 
certainty. But the carrier density value —2.01 x 10^^ cm~^ 
seems to ensure a minimal sensitivity of the Hall resis- 
tance to dissipation. 
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FIG. 5. ARu/Rh as a function of _Rxx ~ Ris 24 x 4^ at rig = 
—2.01 x 10^^ cm^^ (filled blue square), ARn/Rn as a function 
of 7?xx = -R28,34 x ^ at Ms = -2.01 x 10" cm"^ (unfilled 
green square), A7?h/-Rh as a function of 7?xx ~ -Ri8,24 x 
at Tis = —2.1 X 10^^ cm~^ (filled red circle). Errors bars cor- 
respond to measurement uncertainties given within one stan- 
dard deviation, la. 

A similar study was carried out on the ML sample. 
Fig. [5^ shows v = ±2 and i' ~ ±6 Hall plateaus at 
S = 11.7Tandr=1.3K that are typical of half- 
integer QHE in monolayer graphene. Fig. shows 
two couples of Ai?H/-RH = -Rh/(-Rk/2)-1 and i?xx val- 
ues measured with two measurement currents 0.5 fiA 
and 1 fiA at Ug — 6.4 x 10^^ cm^^ on the v — 2 
plateau. Although the deviation strongly increases from 
1=1 /zA, the extrapolation to zero dissipation gives 
Ai?H/i?H(i?xx = 0) = (0.73 ± 2.8) X 10-^. The degree of 
accuracy achieved in the ML and BL samples is therefore 
similar. It is independent of the ratio of the energy gap to 
thermal energy since q 35 kx^b sample is by 2.8 

higher than il"^^^^ in the ML sample. The quantiza- 
tion accuracy is probably determined by the presence of 
the same high concentration of charged impurities in both 
samples leading to the carrier density inhomogeneity and 
low carrier mobility. The impact of charged impurities 
on the QHE breakdown will be discussed in the following 
through detailed analysis of the current dependence of 
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i?xx in the BL sample. 
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FIG. 6. a) Hall (Ra ~ i?06,39) and longitudinal (i?xx = 
Ro6,23 X j^) resistances as a function of ris {v in upper-scale) 
at B = 11.7 T, 7 = 200 nA and T = 1.3 K. b) ARu/Rh as 
a function of 7?xx at ris = —6.4 x 10^^ cm~^. Error bars cor- 
respond to uncertainties given within one standard deviation, 
la. 



IV. DISSIPATION MECHANISM IN THE QHE 
REGIME IN BILAYER GRAPHENE 

A. Current dependence of the longitudinal 
resistance 

Dissipation in GaAs/AlGaAs 2DEG was found to in- 
crease with temperature or current through several mech- 
anisms. At low temperature and low current, carriers can 
backscatter from one edge to the opposite edge through 
localized states by variable range hopping (VRH) with 
soft Coulomb gap, characterized by a temperature behav- 
ior of the conductivity (aovRn/T) exp[-{To{^ioc)/Ty/^] 
where fcBT'o(^ioc) = e^/(47r(roer6oc) and ^loc is the local- 
ization length^""— a lower bound of which is the mag- 
netic length Ib = \/h/eB. Current effect manifests it- 
self as an effective temperature fceT'off = eVnCioc/W^- At 
a higher temperature, conductivity is activated follow- 
ing the behavior ctq exp[— (TAct/r)], where ctq is close 
to e'^/h and weakly dependent on the electron-phonon 
coupling in case of a short-range potential^" but ex- 
pected to be universal and equal to 2e'^/h in case of 
a long-range potential^. Experimentally, VRH mech- 
anism was also observed in monolayer graphene^""— . 
In samples based on exfoliated graphene transferred on 
Si/Si02 substrate, screening of the Coulomb interac- 
tion by the close metallic backgate even restores the 
usual two-dimensional VRH mechanism with a tempera- 
ture dependence exp[— (ri(^ioc)/r)^/^] where Ti(^ioc) ~ 
l/(g(i?F)^ioc) and giEp) is the density of states^^^. Con- 
ductivity activation by temperature was also observed in 
graphene system o^^'^^ . 

On the other hand, there are few reports^iS dealing 
with detailed investigation of the QHE breakdown by 
increasing the current in exfoliated graphene. For semi- 
conductors 2DEGs, several electric field assisted mech- 
anisms have been considered to explain the large in- 
crease of longitudinal conductivity leading to the QHE 



breakdown^: quasi-elastic inter-Landau levels scatter- 
ing (QUILLS) ^^i^° possibly combined with intra Lan- 
dau levels scattering^i^, increase of delocalized elec- 
tron states in Landau levels^^ , ordinary^^ electron heat- 
ing, bootstrap-type^'5' electron heating (particularly ef- 
ficient in large-size samples), and electron percolation be- 
tween sample edges by merging compressible islands§i. 
In a sample made of exfoliated graphene on Si/Si02 sub- 
strate, Singh and co-workers^ deduced from the mea- 
surement of breakdown current dependence on integer 
filling factor that the QHE regime is broken by inter- 
Landau levels scattering in presence of large local electric 
field. 

Fig. [3Ji reports on i?xx dependence on current mea- 
sured at two filling factors v (or Ug values) near ly ~ —4 in 
the BL sample under a magnetic induction of 18.5 T and 
for both temperatures 0.35 K and 1.5 K. It displays ex- 
ponential increases of i?xx over three orders of magnitude 
above a critical current. More precisely, one can define 
a breakdown current 1^ by the value above which con- 
ductivity exceeds 2.10~* S. /c linearly decreases for de- 
creasing V values departing from the filling factor v = —4 
in the range from approximately 1.5 /^A to 0.5 ^A (see 
Fig. [5^). The breakdown current at T = 0.35 K is slightly 
higher than at T = 1.5 K. This behavior is also observed 
in GaAs samples^®. 

Fig. [7] clearly displays the existence of four cur- 
rent regimes for the conductivity calculated by tixx = 
RlJi^lx + Rk)- We will later discuss the first regime 
I for very high currents. For currents down to 1.5 fj,A 
(second current regime II) , conductivity (Txx decreases by 
decreasing the current following a unique fitting function 
(To,,y exp[— Ai?a(j^)/ei?H/] at both temperatures, ctq.^ is 
found quite universal around 0.25e^/h within 30% for 
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v = -4.5 (1=1.5 K) 
v = -4.8 (1=1.5 K) 
v = -4.2 (1=0 35K) 
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FIG. 7. Conductivity (Txx {(? /h) as a function of / ^ at u = 
—4.2 (filled green circle), u = —4.5 (red triangle), i/ = —4.8 
(blue diamond) (filled at T = 1.5 K and unfilled at T = 
0.35 K). Solid lines and dots lines correspond to theoretical 
adjustments at T = 1.5 K and T = 0.35 K respectively. Error 
bars correspond to uncertainties given within one standard 
deviation, la. 
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all V values. Fig. [8^ shows that lS.EJ^v) scales linearly 
with I/, similarly to l^E^-^iv) = (fiwc/2)(l + {v + 4)/2) 
which is the energy difference between the Fermi level 
and the center of the n = —2 Landau level if a constant 
density of states is assumed. The small discrepancy be- 
tween /^Et\i{v) and AEi^{i') results in a deviation of the 
v value for which Ai?a = from —6, center of the n = —2 
LL. The resulting filling factor i-'cdgc can be interpreted 
as the mobility edge which separates localized and ex- 
tended states near the center of the nearest LL. At lower 
current in the third current regime III, <Txx at T = 0.35 K 
decreases more quickly with / and departs from a^x at 
T = 1.5 K that roughly continues to follow the law char- 
acterizing regime IL Reducing the current below 0.7 /iA 
leads to the fourth regime IV where conductivity appar- 
ently saturates at values (Ttm different for the two tem- 
peratures. At T = 0.35 K, the conductivity threshold 
cannot be determined because of the increasing weight 
of some hysteretic charging effect altering measurements 
of (Txx for currents below 0.5 ^A. Consequently, the rea- 
sonable assumption that aT,^ follows the typical tem- 
perature dependence of VRH mechanism cannot be con- 
firmed. 

We remark that the monotonous behavior of 
conductivity following the current dependence 
CTo,i/ exp[— Ai?a(j^)/ei?///] dominates at T = 1.5 K 
in regime II and III over more than three conductivity 
decades down to / = 0.7 /iA (and at T = 0.35 K in 
regime II), and cannot be explained by an activation 
effect caused by a simple heating of electrons by current 
since it does not manifest itself at the lowest temper- 
ature T = 0.35 K in the low current regime III down 
to the same value / = 0.7 /iA. More quantitatively, a 
conductivity increase due to heating by current should 
be described by cto,!/ exp[— A£'th(i^)/fcBT'ci] with Td the 
effective electron temperature resulting from the heating. 
The correct adjustment of data at T = 1.5 K approxi- 
mately above / = 0.7 by ctq,!/ exp[— A£'th(i^)/ei?H-^] 
indicates that k^Td — eVu — eRnI. This would lead to 
an effective temperature Td of 52 K for / = 0.7 fj,A. This 
is not in agreement with the electronic temperature, 
which is obviously close to 1.5 K since conductivity 
starts to be current independent below / = 0.7 /^A, with 
a constant value expected to be determined by the bath 
temperature. Finally, absence of strong asymmetry of 
i?xx values (there are similar within 30%) with respect 
to current direction indicates that there is no strong 
local electron heating in current contact^. This rules 
out any strong role of current contacts 8 (5.9 kfl) and 3 
(< 5.9 kfi) on the observed breakdown mechanism. 

The exponential dependence of conductivity on 
Ai?th(f^) iu current regime II rather directs towards a dis- 
sipation mechanism based on quasi-elastic inter-Landau 
levels scattering (QUILLS) assisted by the electric field. 
In current regime III, the decrease of conductivity at 
T = 0.35 K suggests that the QUILLS mechanism is 
combined with a blockade mechanism manifesting itself 
approximately below T = 1.5 K and a threshold current 



(~ 1 /iA), like a heating mechanism by current. It ap- 
pears that all conductivity curves under regimes II, III 
and IV, at both T = 1.5 K and T = 0.35 K, can be ad- 
justed (see Fig.[7|and Fig.|3Ji) by a unique fitting function 
CTxx = aT,„ + aa^^e:xip[-AEi,{iy)/eRiiI]eyip['~Ec/{kB{T + 
7CTxxVh^))]. cTo^jy ~ 0.25e^/h for all ly values, E^ = 
95 /ieV, 7 = 0.48 K/pW and AE^ near AEth as already 
explained. The contribution (Jt^v is chosen to adjust con- 
ductivity in the low current regime IV only at T ^ 1.5 K. 
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FIG. 8. a) AiSth (meV) as a function of v (red solid line), A_E 
(meV) as a function of v (red filled circle), breakdown cur- 
rents Jc as a function of v a,t T = 0.35 K (blue filled square) 
and T — 1.5 K (blue unfilled square). Error bars correspond 
to uncertainties given within one standard deviation, la b) 
Schematics of inter-Landau levels transitions in case of disor- 
dered 2DEG. 



B. Phenomenological model based on Quasi-elastic 
Inter-Landau Level Scattering (QUILLS) 

In order to explain with the QUILLS mechanism the 
main current dependence of conductivity observed in 
regime II, let us first consider an homogenous electric 
field and harmonic oscillator wave functions for the car- 
riers. The tilting of LLs by the electric field brings 
closer localized states at Fermi energy and extended 
states in the nearest Landau levels, increases the wave- 
function overlap, and thus leads to an increased tran- 
sition probability P between LLs. Through scatter- 
ing processes by phonons and/or charged impurities, 
P is proportional to the wavefunction overlap given 
by exp[~Q^lB^] where Q is the typical direct momen- 
tum between Landau levels Q — ojcB/E (see fig. 8b). 
This should lead to P cx ey.p[-{AEti,{v)/ eElef] = 
ex:p[-{AEth{v)/eVnf{W/lBff-, thus a transition 
probability different from the one observed. But, in 
presence of disorder, at a length scale larger than Z^, 
one expects a dependence exp[— a;/^ioc(j^)] of the local- 
ized state wavefunction tail where ^ioc('^) is the local- 
ization length varying like {y — i/^.)^^^ '^*'^'^-' with the 
filling factor of the Landau level center—. This case 
should be particularly valid in bilayer graphene where 
the energy gap between LLs is very large. From the dis- 
tance X = AEt\i{v)/eE between an initial localized state 
at the Fermi energy and a final extended state in the 
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Landau level at the same energy, one therefore expects 
P oc exp[— Ai?th(i^)/(e-ECioc)]- This model can well de- 
scribe the main exponential current dependence observed 
in current regime II provided that large local electric 
fields with a magnitude around Vh/Cioc are considered. 
Assuming ^loc ^ Ib = 6 nm at = — 4 actually leads 
to a high value of the electric field Vh/^ioc ~ 10^ V/m 
for 1=1 fiA. Besides, the correct adjustment of data 
by this model needs Ai^th to be replaced by Ai?a that 
can be interpreted as the energy difference between the 
Fermi energy and the mobility edge of the nearest LLs. 

We argue that the high concentration of charged impu- 
rities (2.10^^ cm~^) in the substrate can lead to such large 
electric fields. In absence of magnetic field, it was demon- 
strated in section II that charged impurities create carrier 
density fluctuations with a magnitude of 10^^ cm~^ and 
a typical correlation length ^ = 11 nm in the consid- 
ered BL sample notably manifesting themselves as elec- 
tron and hole puddles near the CNP. These fluctuations 
are combined with more macroscopic carrier density vari- 
ations extending over larger spatial scales. Their im- 
pact at high magnetic field in the QHE regime has been 
addressed. Scanning of graphene on Si/SiOj substrate 
by tunneling spectroscopy^^ or single electron transistor 
technique^ has indeed shown that the potential land- 
scape drawn by charged impurities is partially screened 
by the Coulomb interaction and leads to the existence of 
compressible islands surrounded by incompressible strips 
like in AlGaAs/GaAs 2DEG- . Jung and co-authors^^ 
even show that electron/hole puddles at zero magnetic 
field turn into compressible/imcompressible islands in the 
QHE regime. It turns out that the localization length 
^loc, or rather its lower bound /s, as well as the charac- 
teristic length of incompressible strips across which Hall 
potential drops, could be similar to the electron and hole 
puddle correlation length ^. Thus, the existence of large 
local electric field in the BL sample with a typical magni- 
tude Vb_/£, ^ 10^ V/m should result from the strong car- 
rier density fluctuations caused by large concentration of 
charged impurities. Similar explanation was proposed by 
Sing and co-workera-^. Another way to understand the 
impact of the carrier density fluctuations is to consider 
that they turn into spatial variations of the filling factor 
in the QHE regime. Otherwise, the current fiows along a 
path minimizing the dissipation that is expected to occur 
ai V = —4. Given the correlation length of the filling fac- 
tor (or similarly of carrier density) fiuctuations and the 
small width of the sample, it is therefore likely that the 
current flows along a narrow percolating incompressible 
path having a typical width ^ = 11 nm. The potential 
drop concentration across this path leads to the existence 
of large local electric fields. Beyond the enhancement of 
the electric field, the role of charged impurities in the 
QHE breakdown has been investigated in conventional 
semiconductor heterostructure. While charged impuri- 
ties are kept away from the 2DEG by the 10 nm to 40 nm 
thick spacer, acoustic electron-phonon interaction con- 
trols the QHE breakdown because elastic scattering by 



ionized impurities increases the inter-Landau level tran- 
sition rate at higher electric field. But numerical work^^ 
shows that the closer charged impurities are from the 
2DEG the lower the electric field at which they are ef- 
ficient. We therefore propose that a high concentration 
of charged impurities located at only about 1 nm from 
graphene in the BL sample could itself be responsible for 
inter-Landau level transitions, which are in addition en- 
hanced by the strong electric fields introduced by the 
carrier density inhomogeneity these impurities induce. 
This results in low QHE breakdown currents (typically 
0.2 A/m) with respect to expectations in graphene^^ and 
therefore prevents us to observe the VRH regime above 
/ = 0.5 AiA. 

The term exp [-Ec/{kB{T + 7crxxV^H^))] allows the de- 
scription of the temperature effect and the weak heat- 
ing effect by current, clearly visible in current regime HI 
(below 1 /iA). It phenomenologically models a block- 
ade mechanism that can be activated by thermal energy 
above a critical energy Ec = 95 /^eV. The effective tem- 
perature of carriers given hy T* = T + 7crxxVH^ leads 
to the best adjustment of data notably reproducing very 
well the sharpness of the crossover between large and low 
current regimes at T = 0.35 K. Even at T = 1.5 K, this 
exponential term allows a better adjustment of conduc- 
tivity. The proportionality of temperature increase with 
the dissipated electric power means that carriers are very 
badly coupled with phonons of graphene or substrate. 
This results in a large temperature dependence on power 
manifested in the value of the parameter 7 ^ 0.48 K/pW. 
At / = 1 /iA, T* amounts to about 0.5 K, 2 K, and 
5.5 K for V values -4.2, -4.5 and -4.8 respectively. The 
electronic temperature increases all the more so as the v 
value departs from v = —4, because of the higher mean 
conductivity leading to more dissipation. The origin of 
this blockade mechanism manifesting itself at low tem- 
perature and typically clearly visibly below T = 1.5 K in 
our experiment, is not understood. However, it is worth 
mentioning that characteristic energies of phonon absorp- 
tion Hcs/Ib, where Cg is the sound velocity, are 22 meV 
(25.6 K) and 6.6 meV (7.7 K) for phonons of graphene bi- 
layer (c^ = 2 x 10" ms'^) and of Si02 (c^ = 6 x 10^ ms'^) 
respectively, thus well above 95 /leV. On the other hand, 
the energy value 1.8 meV (2.1 K) for phonons of PMMA 
covering the sample (cg = 1.6 x 10"^ ms~^) could be com- 
patible with our observations. An explanation based on 
Coulomb blockade effect in compressible islands is more 
improbable since the low value of Ec would mean over- 
sized islands. 

Finally, the observed disappearance of the exponential 
regime in current regime I (see Fig. |4ti) can naturally 
be explained by the QUILLS mechanism because of the 
overlap integral saturation occurring when Landau lev- 
els are very tilted. At higher currents, resistance slowly 
increases with a polynomial dependence i?xx with 
/? varying between 1/3 and 2/3 for values between -4.8 
and -4.2. 

The linear dependence of Ai?a(i') on demonstrated 
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FIG. 9. a) lni?xx (i?i8,24) as a function of v for four current 
values at T = 0.35 K. Id) A£b = - ln[7?xx/(cro,i. J?ii)]e7?H/ as 
a function of v\ AEa, values are reported as black dot points. 



in the range between v ~ —4.8 and v — —4.2 means 
that density of states is to be quasi constant. It is pos- 
sible to verify this hypothesis on a larger range of fill- 
ing factors i> from the R^x dependence on v measured 
at several currents. Fig. [9^ reports ln(i?xx) as a func- 
tion of v measured with current values 3, 4, 5 /iA for 
which the term exp[— £'c/(fcB(r + 7(TxxV^h^))] ~ 1 has 
no impact. Fig. [^b shows that all four curves displaying 
AEb = -\n[Rxx/{ao,uRii)]eRiiI (with ao,^ = 0.25eyh) 
approximatively merge into a unique curve, except far 
from v — —4. This nicely shows that conductivity well 
follows the current dependence ctq^i/ exp[—AEh{i')/eRHl] 
and reinforces the meaning of Ai?b (i^) as the energy dif- 
ference between the Fermi energy and the mobility edge 
of the nearest Landau levels. A£'b(i^) draws the depen- 
dence of this energy difference on v. This energy reaches 
a maximum value of 45 meV at exactly v = —4, which is 
half the energy gap as expected. Fig. ^jp first shows that 
AEa values deduced at v = —4.2, — 4.5, — 4.8 from 
the adjustment of the current dependence of conductiv- 
ity perfectly match the AEb{v) curve deduced from the 
filling factor dependence of i?xx at different currents. Sec- 
ond, it shows that A£'b(i^) linearity holds on both side 
of = —4 over more than one unit variation. A con- 
stant density of states, as deduced, could be a conse- 
quence of the large Landau level overlap inherent to low 
carrier mobility. Far from v = —4, curves do not su- 
perimpose in a unique curve which means that QUILLS 
is no more the mechanism responsible for conductivity. 
The sub-linearity of ln(i?xx) as observed in Fig. [9^ can 
rather be explained by a saturation of the wavefunction 
overlap at filling factors near mobility edges. Therefore, 
even if we expect an increase of density of states, the 
energy determined near mobility edges in Fig. |9}3 is not 
relevant. Extrapolating the linear behavior of AE\^{i') 
at zero energy should give a reasonable estimate of the 
mobility edge filling factor of the n = —2 Landau level 
^'cdgc(n=-2) = —5.55. This value means that the mobility 
edge energy should depart from the (n ~ —2) Landau 
level energy by 10.35 ± 2.3 meV (120 K). This value can 
be compared with the half width of Landau level pre- 
dicted by the Born approximation to be equal to fi/2Te. 
It matches the lower bound that is calculated equal to 
8.5 meV in considering Tc ~ Ttr = 34 nm in bilayer 
graphene because of the 27r Berry's phase and ignoring 



that Ttr could be larger than Tc because of the long-range 
character of the dominant scattering potential. On the 
other hand, the extrapolation of the linear behavior of 
AE\^{i') at zero energy between i/ = —4 and i/ = leads 
to i'odgo(n=o,i) — ~3, which corresponds to a mobility 
edge shifted from the Landau level center by a larger 
energy of 34.5 meV (400 K). This value cannot be ex- 
plained by the model of the broadening by disorder valid 
for n = —2. 



C. Longitudinal resistance reproducible 
fluctuations 

Fig. Et shows that in the BL sample, at both T = 
0.35 K and T = 1.5 K, the pattern of reproducible fluc- 
tuations of i?xx in the low current regime at / = 0.5 /iA 
is similar to that measured at J = 1 /iA, where con- 
ductivity mainly results from the QUILLS mechanism. 
Measurements at currents above the breakdown current 
2 /iA) have shown a strong decrease of the relative 
amplitude of these fluctuations. From these observations 
we deduce that the QUILLS mechanism adds a conduc- 
tivity contribution that does not itself fluctuate with car- 
rier density. Only the term aT,i/ manifesting itself in the 
low current regime IV has fluctuations. The resistance 
shift due to the QUILLS mechanism by increasing cur- 
rent from 0.5 /zA up to 1 fiA is particularly visible in the 
inset of Fig. |3]; which reports ln(i?xx) as a function of rig 
between —1.9 x 10^^ cm~^ and —2.3 x 10^^ cm^^ away 
from the R^x minimum at rig = —1.8 x 10^^ cm^^. It also 
shows a quite linear relationship between ln(i?xx) and Ug 
for both currents. At / = 1 /iA such a behavior is ex- 
pected since this is a feature of the QUILLS mechanism, 
as also observed at higher current in Fig. [5^1. and Fig.[3]3. 
On the other hand, at / = 0.5 /iA in current regime IV, 
QUILLS cannot account for the linear behavior since it 
is no more the dominant dissipation mechanism, as ob- 
served in current dependence of conductivity in Fig. [T] 
But it turns out that VRH with soft Coulomb gap pre- 
dicts ln(i?xx) (X -(ro(6oc)/r)i/2 cx -Ci^y. Assuming 
^loc oc (i/ — i/c)"^'^, VRH also leads to a near linear be- 
havior of logarithmic conductivity with u, since ln(i?xx) 
should be proportional to (v — i/c)^ '^/^. At = —4.8 
{ris = -2.2 x 10^12 cm-2), VRH would lead to To = 5 K 
and ^loc — 800 nm. VRH can also explain that fluc- 
tuation amplitude decreases as v increases, and that it 
decreases slightly with increasing temperature and more 
strongly with increasing current. This mechanism in- 
deed predicts longitudinal resistance fluctuations result- 
ing from gaussian fluctuations of the localization length 
6oc with an amplitude 51n(i?xx) oc r(-i/2)^(-^3/2) j^^^^ 
that decreases as temperature and ^loc increase. Decreas- 
ing of the amplitude with current can be explained by 
the heating effect by current becoming very significant 
far from v = —4 for currents near 1 //A, as observed in 
regime HI for T = 0.35 K when decreases from —4.2 to 
-4.8 (FigEl): at T*{1 /iA) = 0.5 K at i/ = -4.2 increases 
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up to, T*(l /xA) = 5.5 K at 1/ = -4.2. Thus, the repro- 
ducible fluctuations of i?xx observed are compatible with 
the existence of VRH in the regime IV at low current 
and low temperature. Coulomb blockade in compressible 
islands surrounded by incompressible strips could also 
be considered as a source of conductivity fluctuations. 
In this hypothesis, peaks of conductance would corre- 
spond to the addition of one electron into islands-'^ and 
TrAns?"^ = 1 with Arig the carrier density width of peaks. 
Considering the experimental value of Arig leads to typ- 
ical radius of islands r between 43 nm to 70 nm, thus in 
agreement with values found by others groups^SiZiiZi but 
also not so far from the puddle correlation length at low 
field (11 nm). Although it is difficult to conclude about 
the mechanism at the origin of fluctuations, they can ex- 
plain fluctuations of the Hall resistance i?H observable in 
Fig. |4^ due to the unavoidable residual coupling between 
i?H and i?xx- 

V. CONCLUSION 

To conclude, we have performed quantization tests of 
the QHE in /xm wide Hall bars based on bilayer and 
monolayer exfoliated graphene deposited on Si/SiOj sub- 
strate where electronic transport properties at low mag- 
netic field are mainly governed by the Coulomb inter- 
action of carriers with a high concentration of charged 
impurities. On the Hall plateaus corresponding to Lan- 
dau level filling factor near v — 2 in the ML sample and 
V = —4 in the BL sample, the Hall resistance i?H re- 
spectively agrees with -Rk/2 and -Rk/4 within a relative 
uncertainty of a few parts in 10^, in the limit of zero 
dissipation or at low current in the /iA range. These 
experiments are therefore the most accurate QHE quan- 
tization measurements in exfoliated graphene to date and 
the first in bilayer graphene. At low magnetic field, 
charged impurities mainly located in the silicon sub- 
strate at about 1 nm below the surface and with den- 
sity near 2 x 10^'^ cm^^ reduce mobility, stronger in the 
BL sample < 2300 cm^/V/s) than in the ML sam- 
ple (< 4050 cm^/V/s). These very efficient long-range 
scatterers also induce large spatial fluctuations of car- 



rier density that stays bipolar up to finite density val- 
ues (2 X 10^^ cm~^ in BL). Such density inhomogeneity 
can notably be responsible for the saturation of L$ ob- 
served in the BL sample at low temperature and at fi- 
nite density. In the QHE regime, dissipation leading to 
the QHE breakdown mainly occurs through quasi-elastic 
inter-Landau level scattering (QUILLS) in presence of 
high local electric fields. We claim that a high concen- 
tration of charged impurities very close to graphene ef- 
ficiently assits elastic inter-Landau levels transitions. In 
addition, charged impurities induce a strong filling factor 
spatial inhomogeneity which is favorable to the existence 
of large local electric fields. At low temperature and low 
current, it is observed in the BL sample that dissipation 
also follows an activation law with a typical energy of 
95 /xeV, the origin of which is not understood. As a re- 
sult, breakdown is very anticipated at currents as low as 
1 /iA by enhancement of the inter-Landau level transi- 
tions which prevent from measuring the Hall resistance 
quantization with better accuracy at higher currents. 
This is even more tragic in the small graphene samples 
produced by exfoliation technique. The role of charged 
impurities present in the ML sample is expected to be 
qualitatively the same in the anticipated breakdown, but 
possibly with quantitative differences resulting from par- 
ticularities of the Coulomb potential screening. We then 
conclude that the development of a graphene based quan- 
tum resistance standard able to challenge GaAs would 
require large samples with higher mobility and more ho- 
mogeneous carrier density. To achieve this, the role of 
substrate on which graphene is deposited or grown has 
to be carefully addressed whatever the graphene fabrica- 
tion technique considered. This is a consequence of the 
high sensitivity of graphene electronic transport proper- 
ties to its environment. 
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